Periodic solutions of a resistive model for nonlocal Josephson dynamics 
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A novel method is developed for constructing periodic solutions of a model equation 
describing nonlocal Josephson electrodynamics. This method consists of reducing the 
iy^ . equation to a system of linear ordinary differential equations through a sequence of non- 
^ ^ \ linear transformations. The periodic solutions are then obtained by a standard procedure 
Ch . which are represented in terms of trigonometric functions. It is found that the large 
time asymptotic of the solution exhibits a steady profile which does not depend on initial 
conditions. 
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^ ' 1. Introduction 

The recent studies on Josephson tunnel junctions with high-temperature superconductors 
^ ' reveal that the nonlocal nature of Josephson electrodynamics becomes dominant when 
^ the Josephson penetration depth Aj is shorter than the London penetration depth A^. In 

particular, if we consider a thin layer between two superconductors, the phase difference 
t) across the Josephson junction is described by the following model equation [1-5] 

+ ^j'^v4>t = - sin H f- / Kq — (j)x'x'{x', i)dx' + 7. (1) 

Here, Kq is the modified Bessel function of order zero, uj is the Josephson plasma fre- 
quency, 77 is a positive parameter inversely proportional to the resistance of a unit area 
of the tunneling junction, 7 is a bias current density across the junction normalized by 
the Josephson critical current density and the subscripts t and x' appended to (p denote 
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partial differentiation. When the characteristic space scale Z of is extremely large com- 
pared with Al, the kernel Kq has an approximate expression Ko{x) ~ 7r5{x) where 5{x) 
is Dirac's delta function. Then, the equation (1) reduces to the perturbed sine-Gordon 
equation [6]. In the opposite limit / << A^, one can use the asymptotic of the kernel 
Xodxl) ~ — ln|a:|. In addition, if we restrict our consideration to the overdamped case 
1] » 1 as well as the zero bias current 7 = 0, then unlike the perturbed sine-Gordon 
equation, the equation (1) becomes an integrodifferential (or nonlocal) equation. It can 
be written in an appropriate dimensionless form as 

0t = - sin + H4>,, H<i>, = -P r ^^dx\ (2) 

TT X' -X 

where H is the Hilbert transform operator. Equation (2) may be termed a resistive 
model for nonlocal Josephson electrodynamics [7]. Note that the equation (2) has been 
proposed for the first time in searching intcgrable nonlinear equations with dissipation 
[8]. The general multikink solutions of the equation (2) have been obtained and their 
properties have been investigated in detail [8]. 

In this paper, we report some new results concerning periodic solutions of the equation 
(2). Specifically, we show that the equation (2) can be transformed to a finite-dimensional 
nonlinear dynamical system through a dependent variable transformation. We then lin- 
earize the system of equations to derive a first-order system of linear ordinary differential 
equations (ODEs). Its initial value problem can be solved explicitly to obtain periodic 
solutions. It is shown that the large time asymototic of the periodic solution relaxes to a 
steady profile independent of initial conditions. 

2. Exact method of solution 

2.1. A nonlinear dynamical system 

We seek periodic solution of (2) of the form 

= iln— , f = Y[^sm (3{x - Xj), (3) 

where Xj = Xj(t) are complex functions of t whose imaginary parts are all positive, (3 
is a positive parameter, is an arbitrary positive integer and /* denotes the complex 
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conjugate expression of /. Using a formula for the Hilbert transform, one has = 
— (ln/*/)a;. Substitution of this expression and (3) into (2) gives the following bilinear 
equation for / and /*: 



i{f:f - rft) -\{f- r') - rj - ru (4) 

We divide (4) by /*/, substitute / from (3) and then evaluate the residue at a; = Xj on 
both sides to obtain a system of nonhnear ODEs for xj 

i = 1,2, ...,7V, (5) 

2p n /=i Pi^j - ^i) 

where an overdot denotes differentiation with respect to t. Note that a dynamical system 
corresponding to the multikink solution is derived simply from (5) by taking the limit 
/3 — >■ [8]. As in the multikink case, it will be demonstrated from (5) that the imaginary 
part of Xj remains positive if it is positive at an initial time. 

Before proceeding, it is convenient to introduce some notations: 

V, (6) 
(7) 



(8) 

(9) 
(10) 

Here, Sj,Uj and Vj are elementary symmetric functions of x;,^; and r]i {I — 1,2, N), 
respectively. In terms oi Uj{j — 1,2, ...,N) and Si, / from (3) can be written as 

— Nx — si) 
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Thus, Uj {j — 1,2,...,N) and Si determine the function / completely. In the following 
analysis, we derive a system of equations for Uj. To this end, we find it appropriate to 
rewrite (5) in terms of and r^j as 

where 

Yii^^r^^)-'" = e-^^(-+^^), = no = e''''^. (13) 



^. ^ au^ ^u=i.^^ -2/30, j = 1,2, ...,7V, (12) 



N N 

a = 

Later, it will be shown that a is a constant independent of t and un obeys a single 
nonlinear ODE. 

2.2. Linearization 

Here, we show that the system of nonhnear ODEs (12) can be hnearized in terms of the 
variables Uj defined by (8). We multiply on both sides of (12) and sum up with 
respect to j from 1 to A?^ to obtain 

1 " 

= -^M^ J](-l)Wn-.-2/3tn, 71 = 1, 2, AT, (M) 

^ ^ .=0 

where Ins is defined by 

In deriving (14), we have used the identity 

7„ = 0, -A^+I<n<-1. (16) 
The time evolution of Un follows from (14) with the help of the formulas [9] 

Un = ^ ^(-l)%,tn_„ 1 < n < iV, J^i-iyu.In-j = 0, n > 1, (17) 



^ j=0 j=0 

where uq = 1 and Jq = 1. In fact, differentiating the first formula in (17) by t and 
substituting (14) for in-j, we can show that the quantity /i„ defined by 

hn = Un+-^UNUn —ulf_n + 2PnUn, n=l,2,...,N, (18) 



satisfies the relation 



h„ 



n-l 



n 



n-l 
j=0 



2na 



where 



+n 



(19) 



(20) 



A straightforward calculation using (17) shows that the quantity in the brackets on the 
right-hand side of (20) vanishes identically so that r„ = 0. It follows from this and (19) 
that 

(21) 



h„. = ^ ^(-l)^7i,i„_,-, n = 1, 2, AT. 



n 



j=0 



Note from (13) and (18) that ho = auN/2 — u*^/{2a) = which, combined with (21), 
leads to the relations /i„ = (n = 1, 2, N). Thus, we see that m„ evolves according to 
the following system of ODEs 

-1 



cy ex 

Un + ^«iV«n ^«Ar-n + 2/5nMn = 0, U = 1,2, N. 

It is remarkable that un obeys a single nonlinear ODE of the form 



(22) 



a 2 



+ 2PNun = 0, 



(23) 



and other N — 1 variables ui, U2, un-i constitute a system of linear ODEs. Substituting 
(13) into (23), we can put (23) into a nonhnear ODE for si 



si — —— sinh(2/3 Im si) + iN, 
2ip 



(24) 



where Imsi implies the imaginary part of Si. 



3. Periodic solutions 

The first step for constructing periodic solutions is to integrate (24). It follows from the 
real and imaginary parts of (24) that 
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Resi = 0, Imsi = -— sinh(2/3Imsi) + A^. 



(25) 



Thus, the real part of Si becomes a constant Resi(i) = Resi(O) = b whereas integration 
of the equation for Imsi yields an explicit expression. In terms of a new variable y — 
2/3Imsi, it is given by 

-2/ _ 2^iv (- tanh f + 1) cosh z/Art + {(2/3iV + 1) tanh f - 2/3Ar + 1} sinh z/^t 
^ ~ 2i/jv (tanh f + l) cosh i^^t + {(2/3 AT - 1) tanh f + 2(3N + l} sinh i^^t ' ^ ^ 

where = ^/{PN)'^ + (l/A) and yo = y(0) = 2/3Imsi(0), For n = 1,2, - 1, on the 
other hand, (22) can be written in the form 



u„. = 



(27) 



Note from (13) and Resi = b that a = e~^*^^ becomes a constant. The solution of the 
initial value problem for (27) can be obtained by means of a standard procedure. It can 
be put into the form of a rational function 



«n(t) = %, n=l,2,...,iV-l, 



(28) 



with 



F = 2un (tanh y + l) cosh + |(2/3A^ - 1) tanh y + 2/3 + l} sinh UNt, (29) 
= 2un (tanh y + l) 



Un{0) cosh |/?(A^ - 2n)ii„(0) + ^ti:;^_„(0)| sinh lyj 

(30) 



where = y/ I3'^{N — 2nY + (1/4). We see that the expression (28) with n — N produces 
(26) and hence it can be used for all m„. 

A novel feature of the solution given above will become apparent if one explores the 
large time asymptotic of the solution. Actually, it is easy to see from (28), (29) and (30) 
that as t tends to infinity, Un {n = 1,2, N) approach the following limiting values 



Ur. 



0, n = l,2,...,iV- 1, UN ^ e^'^''{^/A(J3NyT^-2pN). (31) 

The asymptotic form of follows from (3), (11) and (31), giving rise to 

A^7r~ 



~ 2 tan ^ 



V4(/3Ar)2 + 1-1 
2/3 AT 



tern p ( Nx - b 



2/3 



(32) 
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If we introduce a new variable u by u = (p^, then in the hmit t — oo, m behaves hke 



It is remarkable that the asymptotic form of u does not depend on initial conditions 
except for a phase constant b. It represents a train of nonlinear periodic waves with an 
equal amplitude. Since Ui(0) 7^ 0, the initial profile of u has a spatial period tt/P whereas 
that corresponding to (33) is given by n/Nf]. Therefore, as time evolves, there appear 
N identical waves in the space interval ir/jS. The maximum and minimum values of each 
wave are given respectively by Umax = \/^{l3NY + 1 + 1 and Umin = v^4(/?A^)2 + 1 — 1. 
If we define the amplitude of the wave hj A = Umax — Umin, then A = 2, indicating that 
the amplitude becomes a constant independent of the wavenumber. Figure 1 shows a 
typical time evolution of u for N = 2 where the parameters are chosen as /5 = 0.2, Ui(0) = 
g-1.6 e~°'^,M2(0) = e~'^'^, (xi(0) = 4i,X2(0) = 2i). In this example, the wavelength of 
the periodic wave is 7.85. As expected from the asymptotic form (33), we can observe 
two identical waves with an amplitude 2 in the space interval 15.7 at a final stage of the 
time evolution. 



(33) 



^/lijmyTl + (-1)^ cos 2p(Nx - b) ' 



u 



3 



2 



1 
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Figure 1 . Time evolution of u for N = 2 
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4. Conclusion 

We have developed an exact method for constructing periodic solutions of a resistive model 
equation for nonlocal Josephson electrodynamics. Although the basic equation (2) is a 
highly nonlinear equation with a nonlocal term, it can be linearized through a sequence 
of dependent variable transformations. We have obtained various new results. Among 
them, a system of equations (22) is a crucial consequence for the purpose of determining 
the function / given by (11) in closed form. The large time asymptotic of the solution 
exhibits a novel feature. In particular, it relaxes to a steady profile whose functional form 
does not depend on initial conditions except for a phase constant. 

We conclude this paper with some comments. First, we point out that for = 1, the 
periodic solution presented here reproduces an existing solution which has been obtained 
by a different method [10]. The periodic solutions for A^ > 2 appear here for the first time. 
It is interesting to see that the large time asymptotic (32) for satisfies the static version 
of (2), i.e., H(f)x = sin 0. As already pointed out [8], this equation is a model of dislocation 
derived by Peierls [11, 12]. Second, we can perform a similar analysis for a resistive model 
with a bias current 7. It is a relatively easy task to obtain multikink solutions following a 
procedure developed in [8]. Nevertheless, the construction of periodic solutions deserves 
further study. Third, the method of solution developed here is applicable to other nonlocal 
nonlinear evolution equations as well. For example, we will be able to obtain periodic 
solutions of the sine-Hilbert equation H9t — — sin 9, 9 — 9{x, t). It is noteworthy that the 
construction of periodic solutions of the sine-Hilbert equation has been done by a different 
method, but the exphcit solutions have been presented only for A^ = 1 [13]. On the other 
hand, our method will enable us to obtain periodic solutions for general A^. Fourth, the 
periodic solutions can be used to calculate various physical quantities such as the current 
density and the electric and magnetic fields in a Josephson junction. These quantities can 
be compared with experimental results for high-temperature Josephson junctions. The 
solutions to the various problems mentioned above will be reported in a subsequent paper 
as well as a detailed description of the present short communication. 
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